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Abstract. We extend previous studies of mixed-composition extragalactic cosmic-ray source models, by investigating the 
influence of a non-negligible extragalactic magnetic field on the propagated cosmic-ray spectrum and composition. We study 
the transport of charged particles in turbulent fields and the transition from a ballistic to a diffusive propagation regime. We 
introduce a method allowing a fast integration of the particle trajectories, which allows us to calculate extragalactic cosmic-ray 
spectra in the general case, without using either the diffusive or the rectilinear approximation. We find that the main features 
of the mixed-composition models - regarding the interpretation of the ankle and the non-monotonous evolution of the average 
cosmic-ray mass - remain essentially unchanged as long as the magnetic field intensity does not exceed a few nG. 
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1. Introduction 

The transition from galactic to extragalactic cosmic-rays is cur- 
rently one of the most debated issues in high energy cosmic-ray 
physics. On the one hand, the KASCADE experiment (Antoni 
et al., 2005) has shown that the composition is getting heav- 
ier between the energy of the knee (~ 4 10 15 eV) and 10 17 eV. 
The most likely interpretation of this observation is that the 
different elemental components of the galactic cosmic-rays ex- 
perience successive drops (at energies scaling either with the 
charge or the mass) in this energy range. Although this conclu- 
sion seems to be independent of the hadronic model assumed, 
the details of the evolution of the different components remain 
unclear, which prevents, so far, KASCADE observations from 
putting strong constrains on the exact composition of the galac- 
tic component and its evolution in this energy range and above. 
On the other hand, experiments at higher energies (HiRes-Mia: 
Abu-Zayyad et al., 2000; Fly's eye: Bird et al., 1993; Yakutsk: 
Knurenko et al., 2006; Akeno: Dawson et al., 1998; Shinozaki 
et al., 2004; Haverah Park: Ave et al., 2003; HiRes: Abbasi et 
al., 2005) seem to favor a composition getting lighter between 

10 17 eV and 10 18 eV and dominated by light elements above 

10 18 eV. This picture suggests that the transition from a heavy 
galactic component to a light extragalactic component occurs 
somewhere between 10 17 and 10 19 eV. It is however difficult 
to better determine the energy scale and the structure of the 
transition from these experiments (based either on ground ar- 
ray or fluorescence techniques). As a result, several scenarios 
for the transition from galactic cosmic-rays (GCRs) to extra- 
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galactic cosmic-rays (EGCRs) are still viable. In Allard et al. 
(2005, 2007a, 2007b), we proposed to explain the current high 
energy observations with an extragalactic (proton-dominated) 
mixed composition and showed that in this case, the ankle of 
the cosmic-ray spectrum can be interpreted as the signature 
of the transition from GCRs to EGCRs. We also showed that 
the resulting features in the energy evolution of the average 
maximum of longitudinal development of cosmic -ray air show- 
ers, (X max ), are the most compatible with current experimental 
data. These studies were made assuming negligible extragalac- 
tic magnetic fields. 

The extragalactic magnetic fields (EGMF) are currently 
poorly constrained, and often assumed to be negligible. Indeed, 
no convincing mechanism that could create strong magnetic 
fields over very large scales has yet been proposed. However, 
the current upper limits do not exclude the presence of 
magnetic fields of several nG in the extragalactic medium 
(Kronberg, 1994; Blasi et al., 1999). In the case of non negli- 
gible magnetic fields, several effects on the predicted spectrum 
can be expected as a result of the existence of magnetic hori- 
zons (Stanev et al., 2000) evolving with the rigidity of charged 
particles. Strong magnetic fields, for instance, have been in- 
voked to explain the absence of the GZK cut-off (Greisen, 
1966; Zatsetpin and Kuzmin, 1966) claimed by the AGASA 
collaboration (see Aloisio and Berezinsky, 2004; Deligny et al., 
2004). However, a sharp decrease of the cosmic -ray flux above 
~ 5 10 19 eV is observed by HiRes (Thomson et al., 2006) and 
now with larger statistics by Auger (Yamamoto et al., 2007). 

In the context of the second-knee transition model 
(Berezinski et al., 2002; Aloisio et al., 2006), it has been shown 



2 



N. Globus, D. Allard, and E. Parizot: Propagation of UHE cosmic-ray nuclei in extragalactic magnetic fields 



LO 
0.8 


J 4" 


Y^{b 2 ) = 10 nG'. 


LO 
0.8 


0.6 






0.6 


0.4 






0.4 


0.2 






0.2 


0.0 






0.0 



0.0000 0.0001 0.0002 



0.0003 



0.0004 0.0005 






Gauss 


an Fit: 


Amp= 


0.968 +/- 


0.393 


x0= 


0000 +/- 


0.003 




0.006 +1- 


0003 


Cst= 


■0001 +/. 


0393 




v3= 


000 



■0.04 



■0.02 0.00 0.02 

IK || (PO) 



0.04 




-2 

log(r± — r%) (Mpc) 



Fig. 1. Structure of a purely turbulent magnetic field with variance ~<J(\\B 2 \\} = 10 nG and maximum turbulence scale A„ 
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that magnetic fields of the order of a nG could produce a sup- 
pression of the extragalactic component at low energy - say be- 
low ~ 10 18 eV (Lemoine, 2005; Aloisio et al., 2005; Kotera and 
Lemoine, 2007; Berezinsky and Gazizov, 2007) - due to the ex- 
istence of a magnetic horizon caused by the combined effect of 
the magnetic field and the source granularity (see below). The 
importance of this low energy cut mechanism, originally intro- 
duced to limit the predicted extragalactic protons abundance at 
low energy (Lemoine, 2005), is a matter of debate. In the con- 
text of the second-knee transition model, it results in a very 
abrupt change of composition between 10 17 eV and 10 18 eV 
that would induce a very steep (X max ) evolution with energy, 
possibly in conflict with current observations (see discussions 
in Allard et al., 2007, 2007b). Nevertheless, these studies show 
that non-negligible EGMFs could modify the UHECR spec- 
trum. 

In the present paper, we investigate the impact of non- 
negligible EGMFs on the observed cosmic -ray spectrum and 
composition in the case of an extragalactic mixed composition. 
In the next section, we present our modeling of extragalactic 
turbulent magnetic fields and particle transport. We then study 
the different transport regimes of energetic charged particles in 
turbulent magnetic fields. We finally calculate the extragalactic 
spectrum and the energy evolution of (X max ) for different com- 
binations of magnetic fields, source density and source power 
evolution with redshift, and discuss our results. 

2. The modeling of magnetic fields 

Extragalactic magnetic fields are poorly known. Their spatial 
distribution, intensity, coherence length, time evolution and ori- 
gin are uncertain. The presence of fiG fields in the core of 
large galaxy clusters seems to be favored by the observations. 
However, the spatial extension of these large field zones and 
their volume filling factor in the universe are difficult to eval- 
uate. Recent efforts have been made to model local magnetic 
fields using hydrodynamic simulations (e.g. Dolag et al., 2002; 
Miniati and Sigl, 2004; Kang et al., 2007), but their results 
sometimes differ by orders of magnitude. The resulting deflec- 
tions expected for 10 20 eV protons are thus radically different 
under different assumptions and lead to opposite conclusions 
about whether charged particle astronomy is possible or not at 
ultra-high energies. Even though the simulations agree that the 



magnetic field should be very weak in the large inter-cluster 
voids (B S= 10~ 11 G), it is important to note that some po- 
tential sources of magnetic fields in the extragalactic medium 
are omitted in these simulations (see Parizot, 2004; Kotera and 
Lemoine, 2007 for discussions and references therein). In con- 
trast, an interesting simple alternative to complex hydrodynam- 
ical simulations, offering more freedom to test different models 
of magnetic field evolution with the local density, has been pro- 
posed recently by Kotera and Lemoine (2007). 

In view of the above-mentioned uncertainties, we shall use 
here a simplified approach, assuming that the universe is filled 
with a purely turbulent, homogeneous magnetic field, and will 
study the dependance of the high-energy cosmic-ray spectrum 
and composition on several physical parameters such as the 
variance of the turbulent field, the source density and the source 
power evolution. 

To analyze the different propagation regimes of the charged 
particles in a purely turbulent field, we have performed a nu- 
merical simulation of the particle trajectories. The purely turbu- 
lent field is represented by a sum of N m modes as in Giacalone 
and Jokipii (1999) (see Casse et al., 2002; Deligny et al., 2004, 
for an other method based on Fast Fourier Transform): 

N m 

B = Yj Ak ^ n eX P(' k » Z 'n + i3 "), (1) 
1 

where f„ = cos a n x' n + i sin a„y' n , a„ and B n are random phases 
(chosen once for all) and [x' n ,y' n ,z'„] = VR(8 n , <Pn)] X [x,y,z] are 
coordinates obtained by a rotation of the reference frame bring- 
ing the z axis in the direction of the n th contributing wave (i.e. 
k„ is in direction [8„, </>„], also chosen randomly). The ampli- 
tude A{k n ) is determined as a function of ||k„|| according to a 
specific assumption on the type of turbulence (here, we assume 
a Kolmogorov spectrum with wavelengths between A^ m and 
^max) and the chosen field variance (see Giacalone and Jokipii, 
1999, for more details). 

Figure [T] shows some structures of a particular realiza- 
tion of the turbulent field modeled by Eq. ([TJ, assuming a 
Kolmogorov-like turbulence with principal scale A max = 1 Mpc 
and a field variance of V<P 2 II> = 10 nG (hereafter, we refer to 
this quantity as the r.m.s. value, B, or simply as "the magnetic 
field intensity"). Naturally, the values of an individual compo- 
nent of the field are distributed over a Gaussian of mean 0, and 




the coherence length Ac /l ma x/5 — 0.2 Mpc, as expected for a 
Kolmogorov-like turbulence (Harari et al., 2002). 

Once the turbulent field is modeled, the spatial transport of 
charged particles is simply followed by integrating the Lorentz 
equation, = %xB, which we performed with either the 
Bulirsch-Stoer or the forth order Runge-Kutta methods (Press 
et al., 1993 and references therein) with identical results in all 
respects. 

3. Transport of charged particles in turbulent 
magnetic fields 

3.1. Transition from ballistic to diffusive regime 

The transport of charged particles in a turbulent field is concep- 
tually very simple, since, as mentioned above, it comes down 
to integrating trajectories influenced by the sole Lorentz force. 
However, its detailed treatment is limited by the nature of the 
problem, which is essentially chaotic (and of course by our ig- 
norance of the real configuration of the fields). We are thus 
forced to keep to a statistical description of the trajectories. 
To characterize the particle transport and its different regimes, 
we have simulated 100,000 trajectories for protons between 
10 16 eV and 10 20 eV (5,000 trajectories for each energy consid- 
ered) in a 10 nG field (A max = 1 Mpc) (see Fig. [2J. The energy 
loss processes were switched off in these simulations. 

One should first note that individual trajectories only de- 
pend on the particle gyroradius, r g , which can be expressed as 
a function of particle energy E, charge Z and magnetic field B, 
as: 

r g -l.lMpcx|f^. (2) 

In weak fields or at high energy (cf. Eq. |2j, when the particles 
have gyroradii much larger than the coherence length, Ac, they 
are only slightly deflected over that length, with 66(A C ) =s r g /A c . 
Since the direction of the deflections are not correlated from 
one coherence length to the next, particles diffuse in angle (rel- 
ative to the initial direction), with an angular spread <x e ~ t 1 ^ 2 . 
It also appears that the angular decorrelation (i.e, the loss 
of memory of the initial direction) occurs on a longer time- 
scale for higher energy particles, because of their larger rigid- 
ity. Figure [3] shows the evolution of the angular correlation 
£(f) = Ho-w(t) (where «o is the initial velocity vector and «(t) 
the velocity vector after a propagation time t) as a function of 



time, for different proton energies. One can see that the decor- 
relation process spreads over typically one order of magnitude 
in energy, and two orders of magnitude in time. 

Once the particles have lost the memory of their initial 
direction and are completely isotropized, it is appropriate to 
abandon the description in terms of individual trajectories, and 
consider their propagation as a diffusion process, in which the 
average linear distance traveled by CRs increases as Ar 2 = 
6DAt (this equation actually defines the diffusion coefficient, 
D). In the same spirit as in (Casse et al., 2002), we calcu- 
lated the average quantity Ar 2 /6Af - which may be called the 
instantaneous effective diffusion coefficient (IEDC) - for the 
5,000 trajectories at each energy and plot the result as a func- 
tion of elapsed time, Af. Results are shown on Fig. |4] The 
IEDC first increases linearly with time (as CRs propagate in 
straight lines with velocity c), and then a transition occurs to 
a diffusive regime where the IEDC is constant, with a value 
which can be identified with the diffusion coefficient, D(E). 
Of course, the time at which the diffusion regime settles and 
the value of D depend on the CR energy: lower energy parti- 
cles enter the diffusion regime earlier and have a smaller dif- 
fusion coefficient. The diffusion regime can usually be mod- 
eled as a random walk at constant velocity (here v - c) with 
an isotropic redistribution of the direction after a "mean free 
time", which can be identified to the scattering time, r s , after 
which the direction of the particle is (on average) decorrelated 
from the initial direction (see Fig . [3]> . Quantitatively, the diffu- 
sion coefficients are given by the values of the plateau in Fig. [4] 
They are gathered in Fig. [5] showing the smooth transition be- 
tween the low-energy quasi-linear regime, where D(E) oc E 1 ^, 
and the high-energy non-resonant regime, where D(E) oc E 2 , 
in excellent agreement with theoretical expectations. The tran- 
sition between the two regimes of diffusion spreads over an 
order of magnitude in energy, around an energy Eq such that 
2nr g (Eo) = A c . A simple fit of the diffusion coefficient is rep- 
resented on Fig. [5] corresponding to the following formula: 
D(E) = D Bohm (£o)(£/£o) 1/3 + £> B ohm(£i)(£/£i) 2 , where £, 
is such that 2nr„(Ei) = |/l c . The Bohm diffusion coefficient, 
£ ) Bohm(£) = | r g (E)c, is also shown for comparison. 

As shown above, the propagation of CRs in a magnetic 
field smoothly passes from a ballistic regime, where the de- 
flections are small and the distance traveled by a CR away 
from its source grows as r ~ ct, to a diffusive regime where it 
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Fig. 3. Angular correlation, (cos 69), of a set of protons propa- 
gating in a 10 nG turbulent field as a function of time, at differ- 
ent energies indicated by the labels (in eV). 



goes (on average) as r ~ V4Df (we arbitrarily use the numer- 
ical factor 4 for aesthetic reasons). In principle, the diffusion 
regime is always reached, provided one waits long enough (a 
few D/c 2 ). To be specific, one may define Taff(E) = 4D(E)/c 2 
as the time required to reach the diffusion regime. Likewise, 
^diff(£) = cTdiff = V4-DTdiff = 4D(E)/c is the distance traveled 
away from the source before CRs of energy E isotropize and 
diffuse. 

However, CRs may not "survive" long enough to enter the 
diffusion regime. When proton energy losses are included in 
the propagation code, the IEDC curves look like in Fig. [6] A 
diffusion plateau is never reached, because the particles lose 
energy and their instantaneous diffusion coefficient then drops 
to lower and lower values. With the assumed field parameters, 
all particles above 10 20 eV have the same transport properties, 
because they pass below the photo-pion production threshold 
while still in the ballistic regime. 

More precisely, we define a distance, Dgo{E, A), to repre- 
sent the effective "energy-loss horizon" associated with a parti- 
cle that is observed (at the Earth) at energy E and mass A (while 
it was emitted at Eq with mass Aq by its source). Specifically, 
D^)(E,A) is the distance such that 90% of the particles ob- 
served at E and A have travelled a (curvilinear) distance smaller 
than Dtjo from their source. Obviously, this effective distance 
depends in principle on the energy spectrum and cosmic-ray 
composition at the source and on the redshift distribution of 
the sources contributing to the observed flux. 

To estimate Dgo(E,A), we simulated the propagation of 
an extragalactic mixed composition with a spectrum in E" 2A 
(Allard et al., 2005), assuming a uniform source distribution. 
Figure [7] shows the obtained D90 as a function of E/Z for pro- 
tons, He, CNO and Fe nuclei, compared with A^g for magnetic 
fields of 1 nG and 10 nG. The diffusive approximation typi- 
cally holds, for a given value of the magnetic field, whenever 
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Fig. 4. Instantaneous effective diffusion coefficient, Ar 2 /6Af, as 
a function of elapsed time, for protons of various energies, from 
10 16 to 10 20 eV. The magnetic field is turbulent with an r.m.s. 
value of B = 10 nG and principal scale /l max = 1 Mpc. Energy 
losses have been switched off. 



Adis ^ Dlx)(E,A). This depends not only on B and E/Z (as 
would be the case without energy losses), but also on the mass 
of the particle. 

Interestingly, for a 10 nG field, the transition between bal- 
listic and diffusive regimes for protons in the GZK range 
(3 10 19 -10 20 eV), occurs just around the GZK horizon, as il- 
lustrated in Fig. [7] This change of propagation regime in the 
GZK region of the CR spectrum is one of the ingredients of 
a possible alteration of the standard GZK feature, as we now 
briefly discuss (see also Deligny et al., 2004). 

3.2. Magnetic horizons and modification of the GZK 
argument 

To see how magnetic fields can modify the appearance of the 
CR spectrum around the GZK feature, usually calculated with 
B — 0, it is useful to introduce the concept of magnetic horizon, 
i.e. the maximum distance which (most) CRs can travel away 
from their source in a given magnetic field. This is set by the 
energy loss time horizon, estimated here as T90 = Dgo/c - or 
the age of the source, t s , if it is smaller -, and the diffusion coef- 
ficient (or the IEDC if the diffusion regime is not reached). For 
high energy particles, roughly propagating in straight line, the 
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Fig. 5. Diffusion coefficient of protons as a function of energy, 
for a turbulent margnetic field with r.m.s. value B — 10 nG 
and /Imax = 1 Mpc. The red dots are the computation re- 
sults, with asymptotic behaviors indicated by the thin dashed 
lines. The thick blue line is the simple fit proposed in the text. 
The Bohm diffusion coefficient is also shown for comparison 
(dashed line). 

magnetic field has no influence and one simply gets the usual 
GZK horizon. For low energy particles, the diffusion process 
prevents significant propagation beyond the magnetic horizon: 

R lmgn (E) - ^4D(E) min(r 90 (£),fs). (3) 

According to the standard GZK argument, the CR flux above 
10 20 eV should be sharply reduced compared to that be- 
low 3 10 19 eV, say, because of the sudden decrease of the 
GZK horizon distance. However, in the presence of relatively 
strong magnetic fields, low-energy particles also have a limited 
range. In this sense, magnetic fields behave as "low-cut filters", 
whereas the CMB makes the universe a "high-cut filter" (GZK 
effect). The association and tuning of such filters can lead to 
various shapes of the propagated CR spectrum, different from 
the simple "universal" GZK feature. 

In a 10 nG field, for instance, the magnetic horizon should 
severely limit the spatial propagation of low energy particles 
whereas high energy particles (above 5 10 19 eV) propagate in 
straight line. One can then imagine a situation where the same 
number of sources could contribute at 10 19 eV and 10 20 eV! 
This does not mean, however, that there should not be any GZK 
feature in that case. Indeed, while magnetic fields prevent CRs 
from diffusing far from their sources, they also increase the CR 
density around each source, i.e. within the magnetic horizon. 
This is nothing but the magnetic confinement effect, very fa- 
miliar for Galactic CRs. Since the number of particles remains 
the same, both effects exactly compensate, provided that the 
magnetic confinement spheres (of radius ru) centered on the 
different sources are big enough to merge. If this is the case, 
the propagated CR spectrum with and without magnetic field 
are identical: the spatial horizon of the original GZK argument 
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Fig. 6. Same as Fig.W with energy losses switched on. 



simply translates into a time horizon, CRs at 10 eV being 
able to contribute for a much longer time to the observed flux 
than CRs above 10 20 eV. In fact, the propagated spectrum only 
depends on the time-of-flight distribution of the detected CRs, 
and it obviously makes no difference whether the trajectories 
were curved or not (see the formal demonstration of the "prop- 
agation theorem" in Aloisio and Berezinsky, 2004). 

If the magnetic confinement spheres do not fully merge, 
however, interesting spectral effects can appear. The exact 
shape of the UHECR spectrum thus depends crucially on the 
source granularity, i.e. the typical distance between the most 
important sources, A/? s (see Deligny et al., 2004; Aloisio and 
Berezinsky, 2004). The attenuation of the expected GZK fea- 
ture can then be expected for some particular combinations of 
the EGMF and source granularity. Let us note, however, that 
this scenario has been invoked to explain the absence of the 
GZK cut-off claimed by AGASA, but is now less attractive af- 
ter the results of the HiRes (Thomson et al., 2006) and Auger 
(Yamamoto et al., 2007) experiments, finding a highly signifi- 
cant flux suppression above ~ 5 10 19 eV. 

3.3. Effect of magnetic horizons at low energy 

Another effect related to the existence of magnetic horizons 
may arise if we are outside the confinement sphere around a 
source, which is smaller and smaller as the energy decreases. 
In this case, particles emitted by that source simply could not 
reach us. The induced suppression of the flux at low energy 
could be of crucial importance (between 10 17 eV and 10 18 
eV, or even above in the case of strong fields) and impact the 
transition between Galactic and extragalactic CRs. Such a sce- 
nario has been invoked in the context of second-knee transition 
model to lower the contribution of the extragalactic protons 
around 10 17 eV. It has been shown (Lemoine, 2005; Aloisio 
et al., 2005) that, for a typical distance between sources of the 
order of 50 Mpc, a 1 nG field (with /i ma x=l Mpc) could sig- 
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Fig. 7. Effective energy loss horizon, Dgo(E,A), for different 
nuclei as a function of E/Z, compared with the diffusion dis- 
tance A^giE/Z) for 1 nG and 10 nG magnetic fields. 

nificantly suppress the extragalactic proton component below 
10 18 eV and thus avoid a too large proton fraction at 10 17 eV. 

In the next section, we study some potential effects of extra- 
galactic magnetic fields on the cosmic -ray spectrum and com- 
position around the ankle, in the case of an extragalactic mixed 
composition. The influence of the relevant physical parameters 
is examined. 

4. Calculation of the propagated EGCR spectrum 
and composition 

4.1. Building a trajectory database 

Non-negligible magnetic fields make the computation of the 
extragalactic cosmic -ray spectrum conceptually and techni- 
cally more complicated, because there is not anymore a one- 
to-one relation between the distance of a source and the time it 
takes a particle to reach the Earth after it is emitted. When mag- 
netic fields are neglected, one can simply sum the contribution 
of all sources in the universe by integrating over the source dis- 
tance, and for each source distance use the source power at the 
retarded time (i.e. the unique time of emission of the particles 
contributing to the observed flux), allowing for a possible evo- 
lution of the source, and of course apply the relevant energy 
losses, i.e. compute the injection energy of a particle now ob- 
served at energy E. If the extragalactic magnetic fields are not 
negligible, and especially when a diffusion regime can settle, 
high-energy particles emitted at a given cosmic time - i.e. red- 
shift - can be reaching the Earth now, and thus contribute to 
the measured cosmic-ray flux, from sources in a large fraction 
of the universe. 

Note that we make use of the word "redshift" here to repre- 
sent a time, not a distance. Although these notions can be iden- 



tified for a rectilinear propagation (photons or charged cosmic- 
rays with negligible magnetic fields), the most relevant notion 
for the present calculation is time, since it parameterizes both 
the evolution of the sources and the changes in the conditions 
of propagation - most notably the modification of the photon 
environment (CMB and infrared backgroud radiations) which 
control the energy losses of the cosmic rays. We will thus inte- 
grate over the time of emission of the particles, instead of the 
source distance. 

The computation of the EGCR spectrum with magnetic 
fields can be formally simplified if one assumes that the diffu- 
sive regime of propagation holds for all particles at all energies. 
In that case, one can write down the diffusion equation (see 
Aloisio and Berezinsky, 2004 and Berezinsky and Gazizov, 
2007, for a generalization to time-dependent energy losses and 
an expanding universe), and solve it with the appropriate diffu- 
sion coefficients (see above). However, the highest energy par- 
ticles cannot be expected to reach the diffusive regime before 
losing a significant fraction of their energy, even with magnetic 
fields near the highest values allowed by astrophysical obser- 
vations. Using the diffusion equation at low energy and a stan- 
dard rectilinear propagation at high energy would be appropri- 
ate in some cases, but the transition from one energy range to 
the other is likely to contain important information that would 
be missed with such an approach, and furthermore the transi- 
tion energy should be different for different nuclear species. We 
thus chose to avoid approximations (either diffusive or rectilin- 
ear) and compute explicitly the trajectory of the particles in the 
EGMF, whatever their energy and mass, as described in Sect. [3] 

To save computational time, we developed a propagation 
scheme that effectively short-cuts most of the trajectory in- 
tegration. The idea is to build a "database" for the particle 
transport, generated once and for all by integrating many tra- 
jectories with great accuracy over a limited time of propaga- 
tion, and then simulate a realistic trajectory over a much longer 
timescale by putting one after the other the required number of 
(pre-computed) bits of trajectories, drawn randomly from the 
database. Of course, one can in fact forget everything about the 
trajectory itself, since we are only interested in the resulting en- 
ergy, mass and position of a particle (with respect to its source) 
as well as its direction of propagation (with respect to its initial 
direction) after a given time of propagation. Joining random 
bits of trajectories one after the other, one obviously loses any 
correlation between them, so one needs to make sure that the 
bits of trajectory in the database are long enough for such corre- 
lations to be negligible. In practice, this is done by choosing tra- 
jectory lengths larger than the coherence length of the magnetic 
field, A c (see above). We checked that lengths of 1 Mpc and 
10 Mpc gave the same results, and finally chose to divide the 
trajectories into steps of length Z step = 10 Mpc (it takes longer 
to build the database than with 1 Mpc, but then the propagation 
runs to compute EGCR spectra and compositions are quicker). 
Note that we are referring here to curvilinear distances, i.e. the 
trajectories are integrated with full accuracy (using a Bulirsch- 
Stoer integrator to solve the equation of motion, with accura- 
cies from 10~ 8 to 10~ 5 ) over a time f step = l step /c 32.6 Myr. 

The database was built from 10 6 proton trajectories, divided 
into 5,000 trajectories for each of the 200 initial energies con- 
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Fig. 8. Geometry of the trajectories used in the fast integration 
method. The particle is initially injected in O with velocity uo- 
After a time t, it is in M(f), at an angle 6 pos from its original 
direction, and with a velocity u making an angle we \ with the 
radial direction OM (azimuthal angles are chosen randomly as- 
suming isotropy). 



10 2 
10 1 
10° 



\j 10 1 

Q. 



10"' 



10"V • 



10 4 



Bulirsch-Stoer integration 
Fast integration method 



B =10nG 

rms 

X = 1 Mpc 

max 



10 11 



10 1 



10 13 
E(eV) 



10 19 



10 2 



Fig. 9. Comparison between the diffusion coefficients obtained 
with the detailed Bulirsch-Stoer integration (diamonds) and 
with our fast integration method (filled circles). 



sidered, between 10 16 eV and 10 20 eV. The trajectories of nu- 
clei with different charges, Z, are identical if the rigidities are 
the same, i.e. for identical E/Z ratios. Each trajectory is sum- 
marized, after a propagation time t, by a triplet of quantities 
(see Fig. [8}: [r(f), cos <? pos (f), cos vel (f)], where r(f) = \\OM(t)\\ 
is the rectilinear distance traveled away from the source after 
a time f, 6* pos is the global deflection with respect to the initial 
direction of propagation (OM ■ uo = rc cos pos ), and # ve i indi- 
cates the current direction of propagation (with velocity u) with 
respect to the radial vector OM (not with the initial velocity). 
We verified that the corresponding azimuthal angles need not 
be registered for the purpose of our calculations, i.e. we can 
choose them randomly at each step (thus assuming isotropy in 
<£p OS and 9? ve i). 

The distribution of these (triplets of) quantities are tabu- 
lated for each class of trajectories (with a given rigidity), and 
used in the fast integration method as indicated above. An il- 
lustration of the efficiency of the method is given in Fig. [9] 
where we show the diffusion coefficients obtained with the 
fast method, compared with the original ones obtained with an 
accurate integration of the whole trajectories: the differences 
never exceed 15%, which is negligible given our (absence of) 
knowledge of the magnetic fields exact structure and intensity. 

Note that our main goal is to investigate and understand 
the potential effects of EGMFs on the high-energy cosmic -ray 
spectrum and composition, notably at the GCR/EGCR tran- 
sition, not to work out a particular example based on spe- 
cific assumptions about the structure of the magnetic field in 
our region of the universe. Therefore, we simply considered 
a stochastically homogeneous, turbulent magnetic field with 
r.m.s. intensities varying from to 100 nG. 

As a last step to a realistic simulation of the transport of 
high-energy particles in the extragalactic medium, we need to 
plug energy losses and nuclear photo-dissociation into the spa- 
tial trajectories. This is done using the interaction engine pre- 



sented in (Allard, et al., 2006, 2007a, 2007b), by correcting 
the energy of the particle after each time step. This approxi- 
mation was found to be convenient and reliable enough, as the 
potential change of energy during a time step is negligible with 
our parameters. Indeed, the energy (or rigidity) loss time, Ti 0S s, 
of low-energy particles is much larger than the trajectory time 
step, while the highest energy particles with Ti oss > f step have 
essentially rectilinear trajectories and are thus not significantly 
affected by magnetic fields. As for nuclei, energy looses are 
mainly due to photo-disintegration processes, which do not al- 
ter their rigidity over a time f step . We could then check that the 
results of Fig. [6] were reproduced accurately by our fast inte- 
gration method. 

4.2. Propagated spectra from an individual source and 
a source distribution 

The method just described allows us to follow accurately a par- 
ticle in energy space as well as in mass space and geometrical 
space (although only under the assumption of homogeneous, 
turbulent EGMFs), over a large period of time and eventually 
over the age of the universe or of the sources. Finally, to ef- 
fectively compute the propagated EGCR spectrum and com- 
position from a given source model (i.e. a given choice of the 
initial spectrum, initial composition, distribution of sources and 
source evolution), we need to identify which particles actually 
contribute to the cosmic rays observed now and weight the ef- 
fective energy and mass distribution contributed by each source 
according to its power at the time of emission. 

As recalled above, contrary to the case of negligible mag- 
netic fields, distant sources will contribute from a large pe- 
riod of time, not from the sole (rectilinear) retarded time. Each 
source is characterized by its injection spectrum and composi- 
tion: the number of cosmic rays of nuclear species i (with rel- 
ative abundance y, at a given rigidity E/Zj) emitted per second 
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and per unit energy is given by: 



d 2 Ni(E, t) 
dtdE 



= Qi(E, t) = jiQiE/Zt, t) = nQiE/Zi) x /(f), (4) 



where the injection spectrum, Q(E/Z), is assumed to be con- 
stant in time and identical for all species (at fixed rigidity), 
and the dimensionless function /(f) describes the evolution 
of the source power, usually parameterized by the redshift, z' 
f(t) = g(z), where the relation between z and f is given by 
the underlying cosmological model (we assume a standard uni- 
verse with Q m = 0.3, Q A = 0.7 and H = 71 kms^Mpc -1 ). 
The spectrum Q(E) is normalized to an injected power of 
1 erg/s above a minimal energy E^ a : 



Qi(E)EdE = lergs~ 



(5) 



so that the source evolution function f(t) is directly related 
to the source luminosity at redshift z(t), £. s (z), through the 
simple relation /(f) — £, s (z(t))/(l erg s Note that the nu- 
clear abundances at a given energy are related to the "rigidity- 
abundances" y, through the actual shape of the injection spec- 
trum: a, oc yiQ(E/Zj)/Q(E), i.e. a-, oc y,Z: 1 if the EGCR source 
spectrum is a power-law, Q(E) oc E~ x . 

We now define a time-dependent "contribution function", 
F(A,E,r;t), representing the contribution to the observed 
EGCR spectrum of the particles emitted by a source located 
at a distance r, at a given lookback time f (corresponding to 
a redshift z(f)). To compute this function, we use a Monte- 
Carlo procedure, injecting isotropically a large number, No, 
of particles, distributed according to the assumed source spec- 
trum and composition, and propagating them from the injec- 
tion time f to the present time (i.e. from redshift z(t) to z — 0). 
After that time, the particles are spread over a large volume 
around the source, with a variety of masses and energies as 
a result of their interactions with the photon background. We 
then simply register their number density, energy spectrum and 
composition in different bins of distances from the source, 
r, (where isotropy is assumed, according to our simplifying 
assumptions). Specifically, we count the number of particles 
N{Aj, Ej, Tk\ t) found with mass A, in the j' h energy bin, i.e. be- 
tween Ej and Ej+i = Ej + AEj, and in the k ,h distance bin, 
i.e. between r k and r^ + \ = + Ar^, from which we derive the 
contribution function: 



F(Ai, E, r; t) = 



N(Aj,Ej,r k ;t) 
AEj X 4nrlAr k x No ' 



(6) 



for Ej < E < E j + AEj and r k < r < r> + Ar k , where we divided 
by No to obtain a differential density per injected particle. Note 
that F(Aj, E, r; t) has the dimensions of an energy distribution 
function, Eev^cirT 3 , and depends on the injection spectrum 
and composition. 

To obtain a sufficiently precise description of this contribu- 
tion function, we followed no less than 10 10 trajectories (which 
would be impracticable without a fast integration scheme, as 
discussed above), with a variety of nuclei and initial energies 
ranging from 10 17 5 eV to 10 21 eV, and analyzed the particle 
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Fig. 10. Propagated EGCR spectra for pure proton sources with 
an injection spectrum in E~ 2 6 , compared with the HiRes data 
(Bergman et al., 2005), for different realizations of the source 
positions, randomly drawn with a density of n s = 10 5 Mpc~ 3 , 
and a magnetic field of B = 1 nG. The case of a continuous 
source distribution with B — is also shown (dashed line). 



distribution in mass and energy after propagation times rang- 
ing from 10 Myr to 13 Gyr. The procedure was repeated with 
different source spectra and r.m.s. values of the magnetic field. 

Once the "contribution function" is known for a given 
model, the cosmic-ray flux contributed now by an individual 
source at distance R is obtained by a simple time integration 
over the whole period of activity of the source, weighting the 
contribution function at a given look-back time f by the corre- 
sponding number of particles injected per second (since F is 
normalized to one injected particle: see Eq.[6]l. This is equal to 
the source power, X, s {t), divided by the average energy of a par- 
ticle at the source (E) s = £i Qi(E)EdE/ 2/ Qi(E)dE. 
Finally, the contribution of an individual source is given by: 



0(A, E) 



c 







F(A,E,R; f)x =^ x 



df 



<£>s l+z(t) 



(7) 



where the factor c/An simply converts a CR density into a flux 
(in cm _2 s _I sr _1 EeV _I ), and the factor 1 +z(t) takes into account 
the time dilation between the injection time and the observation 
time, associated with the universal expansion. 

Finally, the total EGCR spectrum is obtained by summing 
the contributions of all the sources in the universe. In practice, 
we assume a given source density (number of "standard can- 
dle" sources per Mpc 3 ), draw randomly of set of source dis- 
tances accordingly, and add the contribution of each source 
as given by Eq. ^Jj with the appropriate distance. Note that 
changing the source evolution model does not require to com- 
pute new "contribution functions", since any evolution function 
-£s(z) can be implemented in the last (quick) integration step, 
Eq. 0. 
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Fig. 11. Propagated EGCR spectra for a mixed source compo- 
sition and an injection spectrum in E~ 2A , compared with the 
HiRes data (Bergman et al., 2005), for different source densi- 
ties and EGMF intensities: continuous source distribution (with 
B = 0); n s = 10~ 3 Mpc~ 3 and B - I nG (A max = 1 Mpc); 
n s = 10~ 5 Mpc 3 and B — 1 nG. Each line of a same group cor- 
responds to a different Monte-Carlo realization of the source 
distribution. 



5. Results 

5.1. Pure proton sources 

Figure [10] shows the propagated EGCR spectrum computed 
with the above method in the case of pure proton sources with 
a density of 10~ 5 Mpc -3 , injecting cosmic -rays with an initial 
power-law spectrum in E~ 2 6 in a homogeneous turbulent ex- 
tragalactic magnetic field with B ms = 1 nG and A max = 1 Mpc. 
Twenty different realizations of the local source distribution 
were considered, in order to explore the so-called cosmic vari- 
ance. The propagated spectra are normalized to fit the UHECR 
data from the HiRes experiment (Bergman et al., 2005). We 
also show the results obtained with negligible magnetic fields 
and a totally uniform source distribution (dashed line). As can 
be seen, the extragalactic spectrum is essentially unaffected 
above 10 18 eV, apart from the sharper cut-off at high energy, 
due to the finite distance of the closest source in a realistic 
source model. At lower energy, the expected cosmic -ray flux is 
reduced as a result of the magnetic field horizon effect, and just 
as for the GZK cut-off above 10 19 ' 7 eV, the exact spectrum mod- 
ification is subject to cosmic variance fluctuations, up to a fac- 
tor of 2 in flux. These results are in good agreement with pre- 
vious studies (Lemoine, 2005; Aloisio and Berezinsky, 2005; 
Berezinsky and Gazizov, 2007; Kotera and Lemoine, 2007). 

In the following, we focus on the case when EGCR sources 
inject not only protons, but a mixed composition of nuclei. 



5.2. Mixed composition sources 

As we have shown in previous works (Allard et al., 2005, 
2007a, 2007b), the overall high-energy cosmic ray spectrum 
can be modeled as a superposition of contributions from 



sources distributed throughout the universe, injecting cosmic 
rays with a power-law spectrum in E~ x and a mixed compo- 
sition, similar to what is known from low-energy cosmic rays 
in our Galaxy. The required spectral index depends somewhat 
on the assumed source evolution function. A good fit of the 
data is obtained with a spectral index of x 2.2-2.3 if one as- 
sumes that the source power scales with the star formation rate 
as a function of time. Here, we focus on the influence of the 
magnetic field on the propagated EGCR spectrum, and simply 
assume that the source power remains constant, which favors a 
spectral index of x 2.3-2.4 (we choose x - 2.4, which gives 
better fits in the presence of a magnetic field). 

5.2.1. Magnetic horizons and EGCR contribution at 
the ankle 

Just as in the case of pure proton sources, and for the same 
reasons, the propagated mixed composition spectrum depends 
on the typical source density and on the particular distances 
of the most nearby sources (cosmic variance). This is shown 



on Fig. 1 1 where the EGCR spectra were computed assuming 
an EGMF of B = 1 nG and source densities of 10~ 3 Mpc -3 
or 10~ 5 Mpc~ 3 . The (unrealistic) case of a perfectly uniform 
source distribution is also shown for comparison (in that case, 
we set B — 0, since the magnetic field is irrelevant). As al- 
ready noted in previous studies (Aloisio and Berezinsky, 2004; 
Berezinsky and Gazizov, 2007), it is natural that the propagated 
spectra obtained with increasing source densities converge to- 
wards the "universal" spectrum obtained assuming a contin- 
uous source distribution. For discrete sources, a flux deficit 
is expected at low energy as a result of the energy depen- 
dent magnetic horizon (see above), and for a given magnetic 
field intensity, its amplitude is larger for lower source densities 
(see Lemoine, 2005; Aloisio and Berezinsky, 2005; Kotera and 
Lemoine, 2007 for the study assuming the pure proton case). 



On Fig. 12 we also show the propagated spectra obtained 
for different Monte-Carlo realizations of a source density of 
« s = 10~ 5 Mpc 3 , for different values of the magnetic field in- 
tensity: B = 0, B = 0. 1 nG, B = 1 nG and B = 10 nG. Naturally, 
the magnetic horizon effect (depleting the cosmic -ray flux at 
low energy) is larger for larger field intensities. In the case of a 
10 nG field (close to the upper limit determined by Blasi et al., 
1999, although for a larger value of the field coherence length) 
the EGCR spectrum can be modified up to energies as high as 
~ 10 19 eV for the proton component (and of course more for 
heavier nuclei). Softer source spectra are thus needed to keep a 
good agreement with the data above the ankle, with power-law 
indices around 2.6-2.7 (the HiRes data are shown as an illustra- 
tion). In any case, intense magnetic fields keep the contribution 
of the extragalactic component almost negligible up to 10 18 eV 
and its composition very light above the ankle, since the heavy 
nuclei coming from even the closest sources cannot reach the 
Earth before being photo-disintegrated. Note that this combi- 
nation of values of the source density and the magnetic field 
intensity implies a late transition from GCR to EGCR, i.e, a 
composition dominated by heavy galactic cosmic rays at least 
up to the ankle, which is compatible with classical ankle see- 
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Fig. 12. Propagated EGCR spectra for a mixed source compo- 
sition and different magnetic field intensities: 0.1 nG (x = 2.4), 
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Fig. 13. Propagated EGCR spectra for mixed composition 
sources with density n s = 10~ 5 Mpc -3 , for various spectral in- 
dices and source power evolutions (the magnetic field intensity 
is B = 1 nG). From top to bottom: x = 2.3 and a source power 
proportional to the star formation rate (SFR); x = 2.4 with no 
evolution; x = 2.3 with no evolution. 

narios (see for instance Wibig and Wolfendale, 2004), but not 



favored by current composition analyses (see Sect. 5.2.5 i 



The above results were obtained assuming sources without 
evolution (i.e with a constant luminosity and comoving den- 
sity). Now, a larger source power in the past - as expected, 
for instance, if the source density or intensity roughly scales 
with the star formation rate in the galaxies - results in a larger 
contribution of the sources at lower energy, since the higher 
energy cosmic rays detected now were emitted more recently, 
from less intense sources. Therefore, in general, a harder source 
spectrum (i.e. with less contribution at low energy) is required 
to fit the high-energy CR data if source evolution is taken into 



account. Typically, a source spectral index of 2.2-2.3 instead 
of 2.3-2.4 is found to better reproduce the observed spectrum 
if the source power is proportional to the star formation rate 
(Allard et al., 2007b). 

In a magnetized universe, the source evolution can some- 
how compensate for the magnetic horizon effect. In Fig. 13 we 



show a few EGCR spectra propagated from mixed composition 
sources with a spectral index of x = 2.3, both without source 
evolution and with a source power proportional to the star for- 
mation rate (SFR), as estimated from Bunker et al. (2004). 
Although different at low energy (where the EGCR component 
is subdominant), the latter case and the case of a uniform source 
distribution (i.e. without evolution) with a steeper spectral in- 
dex of x = 2.4 both give acceptable fits of the HiRes data. We 
finally note that increasing the source power or the source den- 
sity (number of sources per comoving volume) does not pro- 
duce the same effect if the magnetic field is non-negligible, 
since a low source density may result in an almost complete 
absence of low-energy EGCRs at Earth, whatever the source 
luminosity, if the closest source lies outside the magnetic hori- 
zon at these energies. In this respect, an evolution of the source 
number density could in principle counterbalance the magnetic 
field suppression effect more efficiently than the evolution of 
the source luminosity. 



5.2.2. Contribution of different source distances 



To better understand the astrophysical origin of the spectral fea- 



tures, we show separately in Fig. 14 the contribution of differ- 
ent intervals of source distances, D, for the above-mentioned 
source densities. The intervals chosen define the following 
"regions" of space around the Earth (note that they do not 
have equal volume nor equal thickness): region 0, for D < 
10 14 Mpc; region 1, for 10 14 Mpc < D < 10 1 8 Mpc; region 2, 
for 10 L8 Mpc < D < 10 2 2 Mpc; region 3, for 10 22 Mpc < 
D < 10 26 Mpc; region 4, for 10 26 Mpc < D < 10 3 Mpc; re- 
gion 5, for 10 3 Mpc < D. As can be seen, the contributions 
to the high-energy flux are very similar for any choice of the 
source density. 

Since the particles propagate almost rectilinearly, the dif- 
ferences in the overall spectrum arise only from the different 
number of very nearby sources, which leads to a more or less 
pronounced GZK feature (De Marco et al., 2006). At low en- 
ergy, however, the spectrum is dominated by the contribution 
of distant sources when the magnetic field is negligible (due to 
the large value of T90), whereas only nearby sources can con- 
tribute significantly if the magnetic field is high enough to settle 
a diffusive propagation regime and confine low-energy cosmic 
rays in the vicinity of the sources. In that case, distant sources 
only contribute at intermediate energies (between 10 18 eV and 
10 19 eV), where the two effects of the magnetic field and of en- 
ergy losses are not strong enough to prevent the particles from 
reaching the Earth (at least for protons). 
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Fig. 14. Contribution of different intervals of source distances (as indicated by the labels and the "distance line") to the overall 
EGCR spectrum, for a mixed source composition and an injection spectrum in E~ 2A . Left: continuous distribution of sources 
(and B = 0). Center: B rms = 1 nG, source density n s = 10~ 3 Mpc~ 3 . Right: Z? rms = 1 nG, source density n s = 1CT 5 Mpc -3 . 
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Fig. 15. Contribution of different intervals of injection redshifts - i.e. actual times-of-fiight - (as indicated by the labels and the 
"injection redshift line") to the overall EGCR spectrum, for a mixed source composition and an injection spectrum in E~ 2A . Left: 
continuous distribution of sources (and B = 0). Center: B lms = 1 nG, source density n s = 1CT 3 Mpc 3 . Right: B lms = 1 nG, source 
density n s = l(T 5 Mpc~ 3 . 



5.2.3. Contribution of different injection times 

Similarly, we show in Fig. [15] the contribution of different 
intervals of times-of-fiight (i.e. redshifts of injection) to the 
overall propagated spectrum. As already stressed, intervals of 
times-of-flight are essentially intervals of distances when the 
magnetic field is negligible. Otherwise, any given source (at 
a given distance) can contribute from a large period of injec- 
tion times, with correspondingly different "instantaneous prop- 
agated spectra". It appears that the relative contribution of the 
different intervals are very similar for all source densities, down 
to low energy. They would actually be totally identical, even in 
presence of a magnetic field, if the source distribution were 
continuous. Nevertheless, in the case of discrete sources, even 
the closest ones cannot contribute at low energy from low in- 
jection redshifts (i.e, for short times-of-fiight), because of the 
increased path length imposed by the magnetic field. This ex- 
plains the observed depletion of the low-redshift contribution 
at low energy, compared to the case without magnetic field (or 
with continuous sources), the effect being of course stronger 
for the lowest source density (or larger values or the magnetic 
field). The depletion at low redshift is directly visible in the 
total spectrum and is stronger when the source luminosity is 



not evolving with time. In the SFR source luminosity evolution 
case, the lower relative contribution of low redshifts to the to- 
tal spectrum at low energy in the negligible field case explains 
the less pronounced depletion in presence of a 1 nG field (see 
above). 

5.2.4. Contribution of the different nuclei 

It is also interesting to study the effect of the EGMF on the 
different types of nuclei contributing to the propagated spec- 
trum. The relevant quantity being the rigidity, oc E/Z, the effect 
produced at a given energy is larger for heavier nuclei. As an il- 
lustration, we study the case of a homogeneous turbulent mag- 
netic field with B rms = 1 nG and a power-law source spectrum, 

£-2.4 



Q(E) 



Figure 



16 1 shows the contribution of different 
mass groups to the overall EGCR spectrum from mixed com- 
position sources, distributed according to a particular Monte- 
Carlo realization of a source density of n s - 10 5 Mpc -3 . The 
magnetic horizon effect is clearly observed at increasing en- 
ergy for increasing nuclei charges. On Fig. 16 1, we show the 



corresponding relative abundances as a function of energy. As 
expected, the extragalactic component consists almost exclu- 
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sively of protons and He nuclei at low energies, where the dif- 
fusion coefficient of heavier nuclei is too small to let them reach 
the Earth in allowed time. As the energy increases, heavier and 
heavier nuclei can come in, resulting in a larger value of the 
average mass per cosmic-ray, (A), the logarithm of which is 
shown on Fig. [To*} :. The continuous increase of In (A) up to 
~ 10 19 eV contrasts with the roughly constant (extragalactic) 
composition expected in this energy range in the case of negli- 
gible magnetic fields. 

At higher energy, EGCRs get lighter as a result of the pro- 
gressive photo-disintegration of intermediate-mass and heavy 
nuclei, starting with the CNO component (He nuclei are photo- 
disintegrated even earlier, while the heavier components come 
in). At the highest energies, EGCRs propagate almost recti- 
linearly and the mass features are independent of the mag- 
netic field. As can be seen, the different energy scales of the 
GZK-like cutoffs of different nuclei cause the composition to 
get heavier over a short energy range, between 5 10 19 eV and 
~ 1.5 10 20 eV (although the upper end of this interval depends 
on the local density of sources). This composition feature is 
similar to that discussed in the case of negligible magnetic 
fields (Allard et al., 2007b). 



900 



The solid line of Fig. 16 shows the value of ln(A) for the 
extragalactic component only. At low energy, a presumably 
heavy galactic component does contribute significantly to the 
observed spectrum (up to at least a few EeV in the mixed- 
composition models), modifying the evolution of the overall 
ln(A) at the GCR/EGCR transition. This is shown by the dashed 
line on the same plot. 

5.2.5. Cosmic-ray composition 

Unfortunately, composition measurements at high energy are 
difficult and the main observable related to composition that is 
currently accessible is the average penetration depth of cosmic- 
ray air showers in the atmosphere. More precisely, when de- 
tecting the fluorescence radiation accompanying the develop- 
ment of CR-induced air showers, high-energy cosmic ray ob- 
servatories can measure the grammage (in g/cm 2 ) along the 
shower axis at which the maximum of development is reached. 
Although there are some fluctuations from shower to shower 
(due to the essential stochasticity of the first interactions in the 
shower development), this grammage - traditionally referred 
to as X nYdx - is on average an increasing function of energy 
and a decreasing function of the cosmic-ray mass (i.e. (X max ) is 
smaller for showers induced by Fe nuclei than by protons). It is 
thus instructing to compute the average value of X max from the 
above propagated compositions and compare its energy evolu- 
tion with the observational data, as proposed in Allard et al. 
(2005b). 

To this aim, an assumption has to be made about the compo- 
sition of the GCR component, which contributes to the overall 
CR composition, and even dominates at low energy. Following 
previous studies and for the sake of simplicity, we assume 
that the GCRs are composed exclusively of Fe nuclei above 
10 17 5 eV, as suggested by the KASCADE results showing a 
cutoff of lighter nuclei between 10 15 eV and 10 17 eV. Although 
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Fig. 17. Average atmospheric depth at the maximum of devel- 
opment of the CR showers, (X max ), as a function of energy. The 
curves show the expectation for each of the spectra displayed 
on Fig. 1 1 (the case of a source density n s = 10~ 3 Mpc~ 



(the case of a source density « s = 11T J Mpc J is 
omitted for clarity), compared with the HiRes-Mia, Fly's Eye 
and HiRes Stereo data. 



some intermediate-mass nuclei (like C, N, O or Si) might still 
be present around 10 17 5 eV, there influence on the results pre- 
sented here should not be dramatic. For instance, the value of 
(I max > is typically 30 g/cm 2 smaller in the case of a pure Fe 
composition than in the case of a pure CNO composition. 
In Fig. 



17 



we show the energy evolution of (X max ) in 
the case of mixed composition sources with a density n s = 
10~ 5 Mpc 3 and a turbulent EGMF with B lms = 1 nG, as com- 
pared to the data from HiRes-Mia (Abu-Zayyad et. al., 2000), 
Fly's Eye (Bird et al., 1993, rescaled by 13 g cirT 2 upward as 
suggested by Sokolsky et al., 2005) and HiRes Stereo (Abbasi 
et al., 2005) experiments. One curve is shown for each EGCR 



spectrum given in Fig. 1 1 as can be seen, the cosmic variance 



is rather small when one considers an average quantity such as 
(I max ). The case of a continuous source distribution without 
magnetic field is also plotted as a reference (dashed line). 

The actual value of X m - dx for a simulated cosmic-ray shower 
depends somewhat on the underlying hadronic model assumed 
to compute the high-energy cascade. However, the energy de- 
pendence of X max is similar for all models, so differences come 



down essentially to a global shift of the curves in Fig. 17 typi- 
cally by ~ 20 g cm~ 2 (see discussions in Allard et al., 2007a,b). 
Here, we used the QGSJet-II hadronic model (Ostapchenko, 
2004), which leads to the best agreement of the mixed compo- 
sition model with the data. 

As can be seen from the comparison with the no-field case, 
the presence of a 1 nG extragalactic magnetic field does not 
modify significantly the shape of the (X max ) evolution with 
energy. In particular, they leave unchanged the characteristic 
features of the mixed composition model to which the data 
lend support (see Allard et al., 2007a,b). Nevertheless, mag- 
netic fields cause the (X max ) evolution to be slightly steeper 
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log 10 E (bV) log 10 E (eV) log 10 E (eV) 

Fig. 16. Left: contribution of different nuclei to the overall EGCR spectrum, for a mixed source composition (see text), an 
injection spectrum in E~ 2A , a homogeneous turbulent magnetic field with B rms = 1 nG, and (one particular realization of) 
a source density n s = 10~ 5 Mpc~ 3 . Center: corresponding relative abundances as a function of energy (ranges of nuclei are 
identified through their charge number). Right: logarithm of the corresponding average mass, ln(A), as a function of energy 
(solid line). Dashed line: resulting value of ln(A) when a pure Fe nuclei Galactic component is added. 



in the Galactic -to-extragalactic transition region (i.e. between 
jq17.5 e y an( j g_4 jq18 e y^ as a jjjjgpf consequence of the 

EGCR flux reduction (magnetic horizon effect). The effect on 
C^max) remains rather small, however, since the EGCR com- 
ponent does not contribute much at low energy (say, below 
10 18 eV), even with 5 = 0. 

The magnetic suppression of the heavy nuclei in the EGCR 
component also causes (X max ) to be larger when the extragalac- 
tic component dominates, i.e. above ~ 3 10 18 eV, until the bal- 
listic regime of propagation replaces the diffusive regime (i.e. 
above ~ 10 19 eV in this case, for B = 1 nG). As heavier nu- 
clei come in, overcoming the magnetic suppression, the EGCR 
composition gets heavier progressively (see Fig. [To*}:), which 
results in a flatter evolution (smaller increase) of (X ma!l ) than in 
the case without magnetic field. This makes the subsequent in- 
flection point even more pronounced, when the heavy compo- 
nents start falling off (because of the photo-disintegration suf- 
fered by nuclei in the extragalactic photon background) and a 
second transition from heavy to light cosmic rays takes place. 
Note that the second bump in Fig. 16: also produces a vis- 
ible feature in the (X max ) evolution plot, although with the 
same amplitude as in the case without magnetic field, since the 
propagation is essentially rectilinear at these energies. Finally, 
the slight difference at the highest energies is an artefact of 
the continuous source distribution assumed for the reference 
curve (dashed line). In this energy range, the measured spec- 
trum and composition is directly dependent on the few clos- 
est sources in the Earth vicinity, and should be associated with 
strong anisotropies which leave the study of the overall spec- 
trum (gathering the contribution of all sources) without much 
purpose. 

In sum, we find that the characteristic features of a mixed- 
composition model - including a smooth transition from heavy 
to light cosmic rays across the ankle, an inflection point in the 
(X max ) evolution around 10 19 eV (and a small second bump 
below 10 20 eV) - are not attenuated by the extragalactic mag- 
netic fields, but rather accentuated. We furthermore checked 
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Fig. 18. Logarithm of the CR average mass, ln(A), as a func- 
tion of energy, for different values of the magnetic field (as 
indicated), a source spectrum in E~ 2A and a source density 
n s = l(T 5 Mpc- 3 . 



that the features we predict would be even more accentuated by 
increasing the heavy nuclei relative abundances at the sources 
without affecting the agreement of the predicted spectrum with 
data. 

This is a consequence of the dependence of the average CR 



mass on the intensity of the magnetic field, as shown in Fig. 18 
where we plot ln(A) as a function of energy for three different 
values of the magnetic field: 0.1 nG, 1 nG and 10 nG. In the 
latter case, the suppression of the low-energy extragalactic nu- 



clei is so strong (see Fig. 12 1 that the resulting CR composi- 
tion sharply drops around the ankle from the Galactic Fe nuclei 
to an almost pure proton extragalactic component, with only a 
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Fig. 19. Same as Fig. 17 for a 10 nG magnetic field. 
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Fig. 20. Same as Fig. 17 for the spectra displayed on Fig. 
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with a source power evolution proportionally to the star forma- 
tion rate. 



small fraction of nuclei coming in at higher energy. The result- 
ing evolution of (X max ) with energy then appears to be much 



too steep in this energy range, as illustrated in Fig. 19 



6. Discussion 

Following on the study of the spectral and composition fea- 
tures associated with mixed-composition EGCR source models 
(Allard et al., 2005; 2007a; 2007b), we investigated here the 
influence of a non-negligible extragalactic magnetic field and 
found that the main features of such models remain essentially 
unchanged as long as the magnetic field intensity is smaller 
than a few nG. The astrophysical interpretation of the ankle as 
a transition from the Galactic to the extragalactic cosmic-ray 
component remains identical, and relatively strong magnetic 
fields can even require the GCR component to be extended up 
to higher energies, into the ankle. Likewise, extragalactic mag- 
netic fields do not alter the characteristic feature in the energy 
evolution of X max (the main composition observable in cur- 
rent experiments), namely the flattening of the X m . dx evolution 
and the existence of an inflection point between 10 18 5 eV and 
10 19 2 eV, in good agreement with the observations. In most 
cases, it is actually enhanced. For EGMF intensities larger than 
a few nG, the resulting energy evolution of X m - dx becomes in- 
compatible with the data. 

The magnetic field model used in the above calculations - 
namely, a homogeneous turbulence - is admittedly simplistic. 
It is likely that the EGMF shows various levels of structures 
related to the structures in the matter distribution. However, 
in the case of external sources, such structures would only 
lead to smaller propagation effects, since the filling factor of 
the larger fields would be smaller. Note also that Kotera and 
Lemoine (2007) considered more realistic margnetic field con- 
figurations, with intensities between 0.1 and 10 nG, and found 
only a weak dependence of the EGCR (pure proton) spectrum 
above 10 17 eV on the spatial distribution of the field (as con- 
trolled by the specific model relating the magnetic field to the 
matter density). A different conclusion could be reached, how- 
ever, if most sources were located inside highly magnetized 
regions. In this case, a larger suppression of the low-energy 
EGCRs could be expected, as recently discussed by Sigl and 
Armengaud (2005) and Sigl (2007). It should be noted, how- 
ever, that the influence of a magnetic confinement around the 
sources on the propagated spectrum remains uncertain, since it 
greatly depends on the simulation used to compute the mag- 
netic field intensity and its spatial distribution. For instance, 
the magnetic fields simulated in Dolag et al. (2004) are likely 
to have a much weaker effect on the charged particles propaga- 
tion than those of Sigl et al. (2004). 

In the absence of any clear prediction regarding their time 
evolution, we have also assumed constant extragalactic mag- 
netic fields. This should not have a strong impact on our re- 



As noted in Sect. 5.2.1 the effect of a magnetic field at suits, however, since we showed on Fig. 15 that most of the 



low energy can be roughly compensated by an increase of the 
source power at high redshifts (which most contribute to the 
low energy flux). The resulting effect on the energy evolution 
of (^max) is illustrated in Fig. 20 where larger values of X max 
(i.e. a larger fraction of protons, from the EGCR component) is 
observed at low energy, in better agreement with the data and 
essentially identical to the case with a steeper source spectrum 



andfi = 0(cf. Fig. 17 1 



particles contributing to the spectrum above 10 175 eV were 
emitted at redshifts smaller than 0.5, i.e. over a period where 
such an evolution should be quite limited. Even if stronger 
magnetic fields were present at higher redshifts, the results pre- 
sented here would remain essentially unchanged. Note that the 
time evolution model used by Gazizov and Berezinsky (2007) 
also implies smaller coherence lengths at large redshifts, which 
partly compensates the effect of a stronger field. 
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We have also assumed that all EGCR sources are essentially 
identical, i.e with the same injection spectrum and intrinsic lu- 
minosity. This need not be the case, and one may object that 
sources contributing at the highest energies might be less nu- 
merous than lower-energy sources. In this case, the magnetic 
suppression of low-energy particles would be reduced, and the 
results closer to those obtained earlier in the case of a negligible 
magnetic field. 

Finally, we stress that the measurement of the UHECR 
composition is an important step towards the clarification of 
the source phenomenology as well as the progressive transi- 
tion from a Galactic to an extragalactic cosmic-ray component. 
The preliminary Auger data (Unger et al., 2007) appear in good 
agreement with the Fly's Eye results above 4 10 17 eV, and thus 
with the X m - dx features associated with the mixed composition 
models, including a possible flattening of the (X max ) evolution 
at the highest energies. It is however too early to determine 
whether this feature (if any) is associated with the GZK cut-off 
in the spectrum and whether the composition is getting lighter 
again, down to protons, above 10 20 eV. Future Auger data, with 
enhanced statistics and an extension at lower energy (Medina et 
al., 2006), should be decisive in allowing one to distinguish be- 
tween the different GCR/EGCR transition models. In addition, 
as the recent results from KASCADE (KASCADE collabora- 
tion, 2007), KASCADE-Grande (KASCADE-Grande collabo- 
ration, 2007) and Auger (Engel et al., 2007) have shown, a pre- 
cise measurement of the properties of the cosmic -ray air show- 
ers can help better constrain hadronic models at high-energy 
and give more confidence in future composition analyses. And 
since extragalactic magnetic fields may also modify the energy 
evolution of the CR composition, as was shown above, let us 
finally note that new constraints on the strength and structure 
of the EGMF should also result from the analysis of UHECR 
arrival directions and large statistics anisotropy studies, as ex- 
pected from the Pierre Auger Observatory and its Northern 
hemisphere extension (Nitz et al., 2007). 
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